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ABSTRACT 



Context. The measurement of Doppler velocities in spectroscopic solar observations requires a reference for the local frame of rest. 
The rotational and radial velocities of the Earth and the rotation of the Sun introduce velocity offsets in the observations. Normally, 
good references for velocities are missing (e.g. telluric lines), especially in filter-based spectropolarimetric observations. 
Aims. We determine an absolute reference for line-of-sight velocities measured from solar observations for any heliocentric angle, 
calibrating the convective line shift of spatially-averaged profiles on quiet sun from a 3D hydrodynamical simulation. This method 
works whenever there is quiet sun in the field-of-view, and it has the advantage of being relatively insensitive to uncertainties in the 
atomic data. 

Methods. We carry out radiative transfer computations in UTE for selected Ci and Fei lines, whereas the Can infrared lines are 
synthesized in non-LTE. Radiative transfer calculations are done with a modified version of Mum, using the snapshots of a non- 
magnetic 3D hydrodynamical simulation of the photosphere. 

Results. The resulting synthetic profiles show the expected C-shaped bisector at disk center The degree of asymmetry and the line 
shifts, however, show a clear dependence on the heliocentric angle and the properties of the lines. The profiles at /j = 1 are compared 
with observed profiles to prove their reliability, and they are tested against errors induced by the LTE calculations, inaccuracies in the 
atomic data and the 3D simulation. 

Conclusions. Theoretical quiet-sun profiles of lines commonly used by solar observers are provided to the community. Those can 
be used as absolute references for line-of-sight velocities. The limb effect is produced by the projection of the 3D atmosphere along 
the line of sight. Non-LTE effects on Fe i lines are found to have a small impact on the convective shifts of the lines, reinforcing the 
usability of the UTE approximation in this case. We estimate the precision of the disk-center line shifts to be approximately 50 m s"', 
but the off-center profiles remain to be tested against observations. 

Key words. Convection - Sun: photosphere - Line: fomation - Radiative transfer - Sun: granulation - Hydrodynamics 



1. Introduction 

The need for high spatial and spectral resolution in solar obser- 
vations stimulates new advances in telescope and instrumenta- 
tion technology, adaptive optics and image reconstruction tech- 
niques. However, the accuracy of the measurements is often lim- 
ited by calibration issues. This is evident when using Doppler 
shifts of spectral lines to acquire line-of-sight velocities in the 
solar plasma, where the task of finding a local standard of rest is 
often problematic. 

The current work is partly motivated by the install ation of 



the C Risp Imaging Spectropolarimeter (CRISP, Schar mer et al 



2008 ) at the Swedish 1-m Solar Telescope (SST. IScharmer et al 



20031) . Instruments like this usually lack laboratory wavelength 
references, while the use of telluric lines has limitations as dis- 
cussed below. The obvious (and commonly used) solution of let- 
ting the spectral line of interest - averaged over a large enough 
area and long enough time - define a local frame of rest is con- 
founded by the fact that convective motions leave strong finger- 
prints on any spectral line formed in the solar photosphere. The 
statistical average of bright blueshifted and dark redshifted pro- 
files formed in granules and intergranular lanes respectively, pro- 
duces the bl ueshifted C-shap ed bisectors of photospheric lines as 
reviewed bv lDravinsl(ll982l) . 

Send offprint requests to: J.d.l.C.R.: e-mail: jaimeOastro .su.se 



Below, we discuss some methods that have been used to de- 
fine a velocity reference for solar observations. 

- Sunspot umbrae are common standards (e.g.. iBeckersll 1 9771 : 
ISc harmer et"ani2008t lOrtiz et al.ir2010l) . but this only works 
when a suitable sunspot is in the field-of-view and relies 
on the assumption that the umbra is at rest because of the 
convective motions being suppressed by the strong magnetic 
fields. Furthermore, the line must be measurable in the um- 
bra and wavelength shifts which are induced by blends from 
lines uniquely formed in the umbra (e.g. molecules), must 
be accounted for. Finally, observations of umbrae are always 
plagued by stray light from the much brighter quiet photo- 
sphere. 

- Telluric lines can be used to define a laboratory-frame 
reference, which c a n the n be converted to th e Sun . 
iMartinez Pillet et all (Il997h and iBellot Rubio et alJ (l2008h 
used telluric lines to derive the wavelength scale for their 
observations, then converted to an absolute wavelength scale 
on the Sun given the ephemeris constants, time of the obser- 
vations, solar rotation and the laboratory wavelengths for all 
observed spectral lines. 

- A spectral atlas, most commonly the atlas acquired 
with the Fourier Transform Spectrometer at the McMath- 
Pierce Telescope (her eafter called the FTS atlas) of 
iBrault & Neckel (Il987h . can be used to calibrate observa- 
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tions. This was done by iLangangen et al.l (l2007h . who used 
the laboratory wavelength to define the reference for ve- 
locities, assuming an absence of systematic errors in the 
FTS atlas itself and correcting for the gravitational redshift 
(633 m s This method only works at disk center where 
the atlas was recorded. 
- Numerical models can be used to predict convective line 
shifts and thus provide a velocity reference. A two- 
compon ent model of the solar photosp here has been de- 
rived by iBorrero & Bellot Rubior(l2002li from the mversion 
of Fei spectral lines. Since it does not contain any hori- 
zontal velocities, it then is in principle limited to calibrat- 
ing disk-center observ atio ns, for which it has been used b y 
iTritschler et al.l (|2004|) andlFranz & Schlichenmaied (l2009l) . 
However, iBellot Rubio et al. used it with the empiri- 

cal results of Balthasai (119881) to also estimate Imeshifts off 
solar center ILangangen et al. (l2007h employ a 3D numeri- 
cal granulation model to compute the convective blueshift of 
the C 1 5380 A line and calibrate their observations. This ap- 
proach was necessary because the laboratory wavelength of 
the C I line is not known with enough precision to use the 
atlas calibration that is mentioned above. 

We note that telluric hnes cannot always be observed with 
the solar diagnostic lines and that when this is the case, ob- 
serving them can be expensive. Observations with tunable fil- 
ters require high cadence because of seeing variability and so- 
lar evolution. Any time spent recording telluric lines thus de- 
creases the spatial resolution and the signal-to-noise of the sci- 
ence data. This would also apply to any laboratory spectral-line 
source used for calibration. Furthermore, approaches using tel- 
luric lines (or other calibration lines) and spectral atlases all re- 
quire very accurate atomic data to define the zero velocity point 
and thus convert the wavelength scale into velocity scale. Such 
data are unfortunately not available for many of the lines that 
ar e extensively used in sol ar physics as is apparent in the work 
of ILangangen et akl (|2007|) described above. Their approach is 
proposed in this study: using realistic 3D numerical simulations 
of the solar photosphere to compute spatially-averaged line pro- 
files. As the radiative transfer computation is carried out in the 
local frame of rest of the Sun and the adopted atomic data is 
known, the velocity shifts of the spatially-averaged profile can 
be computed accurately for different heliocentric angles, relative 
to the assumed center of the line. This method relies on the de- 
gree of realism of the 3D simulation and requires a large statisti- 
cal sample of spectra to compute the spatially-averaged profile. 
The aim of this work is to provide the community with a set 
of theoretical spatially-averaged quiet-Sun profiles of lines com- 
monly used in solar physics, that are computed for a range of 
heliocentric angles. 

In Sect. |2] the atomic data and radiative transfer details are 
discussed. The 3D HD simulation is described in Sect. |3] In 
Sect. |4] the results are presented together with an error analy- 
sis and a real filter-based observation is calibrated using the new 
data. Section |5]describes the electronic data made available and 
Sect.|6]summarizes the main conclusions of this work. 

2. Line list and radiative transfer 

The calculations were carried out with a selection of lines 
(Table [TJ that are extensively used in solar physics, including 
lines that can currently be observed with CRISP at th e SST. A 
modified version of the radiative transfer code Mum (ICarlssonI 
Il986h was used to compute synthetic spectra from the simula- 



Table 1. 

length. 



Atomic data, sorted by element in increasing wave- 



^0 (A) 


E, 


^og gf 




a<2) 


log 6<3) 


C I 5380.3370 


7.6850 


-1.824 


1238 


0.229 


8.53 


Fe I 5250.2089 


0.1210 


-4.938 


207 


0.253 


7.44 


Fe I 5250.6460 


2.1980 


-2.181 


344 


0.268 


7.71 


Fe I 5576.0888 


3.4300 


-1.000 


854 


0.232 


7.69 


Fe I 6082.7106 


2.2230 


-3.573 


306 


0.271 


7.45 


Fei 6301.5012 


3.6540 


-0.718 


834 


0.243 


7.55 


Fe I 6302.4940 


3.6860 


-1.203 


850 


0.239 


7.50 


Fei 7090.3835 


4.2300 


-1.210 


934 


0.248 


7.63 


Can 8498.0134 


1.6920 


-0.6391 


291 


0.275 


6.31 


Can 8542.0532 


1.7000 


-0.4629 


291 


0.275 


6.31 


Can 8662.1605 


1.6920 


-0.7231 


291 


0.275 


6.31 



Notes. (1) the cross-section in units of the Bohr radius, (2) dimension- 
less parameter used with (1) in the computation of coUisional broaden- 
ing tBarklem et al...l998h . (3) the fitted value of the abundance for each 
line. 



tion snapshots. The Fe i and C i lines were calculated assuming 
local thermodynamical equilibrium (LTE), which is a reason- 
able assumption for the wings of most p hotospheric Fe i lines 
(iGadun et all [19991: lAsphmd et al.ll2000h . iFabbian et aP (l2006l) 
found that non-LTE effects for C i lines are expected to be small 
with the effect on derived abundance being < 0.05 dex. Possible 
errors produced by non-LTE effects will be further discussed in 
SectiongXH 

The Ca ii lines were computed in 1 .5D non-LTE - where the 
n on-LTE problem is solv ed in ID - using the same model atom 
as lLeenaarts et al.l (l2009l) . which consists of 5 bound levels plus 
a continuum. 

Atomic data were obtained from the VALD-2 database 

dKup ka et al. 2000). T he treatment o f van der Waals broadening 
follows that of lBarklem et alj ( 1998h and the required cross sec- 
tions were taken from Barklem et alj (|2000|) . except for the C i 
line provided by Paul Barklem (private comunication). The ef- 
fective quantum numbers «* for this transition are 1.95 3.27. 
The treatment becomes less reliab le when m > 3 .0 but it is more 
accurate than the commonly used lUnsoIdl (Il955h formulation. 

To achieve the best possible fit between the synthetic line 
profile at yU = 1 and the FTS atlas, the elemental abundance 
is individually adjusted for each line. The resulting abundance 
values (last column in Table [TJ should not be regarded as such. 
Instead, they act as free parameters that compensate for errors 
from atomic data, the LTE approximation, and the models them- 
selves. The hope is that by fitting the profile the correct layer of 
the model with the appropriate velocity field is being sampled. 
Since non-LTE effects could be expected to affect the cores of 
strong lines, only the wings of each line were used for the fit of 
the abundance; however, the results change marginally when the 
full profile is used to fit. 

The results for the Can infrared triplet lines must be ap- 
proached with extra caution since the lines are very strong and 
the 3D simulation does not extend all the way up to the chro- 
mosphere. Therefore the calculated core of the line cannot be 
compared to real observations, but the photospheric origin of the 
wings sugge sts the usability of the models with these lines, as 
discussed bv lShine & Linskvl(ll974h . 

We have included blends affecting the Fez and Ci lines. 
The atomic parameters for these were taken from the VALD-2 
database (Kupkaetal. 2000). The blends are weak, but have a 
small effect on the upper part of some bisectors. The Can lines 
have some strong blends in the wings which greatly affect the 
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Fig. 1. Mean bisector velocity as a function of the number of 
snapshots used to compute the spatially-averaged spectrum. 



bisectors, so these blends were not present in the present Ca ii 
calculations. 

3. 3D hydrodynamical photospheric models 

The 3D hydrodynamical simulation consists of 90 snapshots and 
cover s around 45 minute s of solar time (Trampedach et al., in 
prep.: |Pe"reira et al1l2009l) . The simulation does not include mag- 
netic fields. The radiative energy transport was computed with 
an improved multi-group op acity binning of 12 bin s, when com- 
pared with the 3D models of lAsplund et al.l (|2000|) . The original 
simulation of 240 x 240 x 240 grid-points covering a region of 
6x6x4 Mm was interpolated to a finer vertical depth grid in 
the radiative zone and sampled to a coarser grid in the horizon- 
tal plane for radiative transfer calculations. The resulting snap- 
shots had a size of 50 x 50 x 82 grid points representing a 6 Mm 
square. Our tests show that the effect of such a resampling in the 
horizontal plane has a negligible effect on the spatially-averaged 
profiles of the lines. 

The presence of oscillations in the simulation broadens the 
time-averaged profile compared to the individual snapshots. The 
oscillations have a quasi-periodic behavior, so at least one period 
should be included in the time average in order to avoid system- 
atic shifts of the average profiles. Figure [1] shows the average 
velocity (integrated along the bisector) of the time averaged pro- 
file as a function of the number of snapshots used in the time 
series. The limited sample of granules and intergranular lanes 
enclosed in a 6 Mm box reinforces the need to use more than 
one period in order to converge to a stable average for the line 
shift. Therefore, the whole time series was used for computing 
the average line profiles. 

To produce profiles off disk center, the properties of the 3D 
snapshots are slanted assuming periodic boundary conditions, 
and the velocity vector is projected into the new line-of-sight. 
The height scale is also modified to match the new geometry of 
the 3D model. 

4. Results 

In Fig. |2] the resulting line profiles for all heliocentric angles 
are plotted together, along with their bisectors. The C i 5380 A 
line shows the strongest convective blueshift, an expected result 
given that the line is formed very deep due to the hight excitation 
potential of the transition. In deep layers, vertical velocities are 



high and the surface of the intergranular lanes becomes smaller. 
The blueshift decreases monotonously towards the limb. 

A variety of bisector shapes and fi dependences are found 
in the Fei lines. At /i = 1, the lines have the expected 
blueshifted convective C-shaped bisectors. When going off cen- 
ter, the blueshift for many lines increases, while the blueshift 
decreases again at even higher heliocentric angles, sometimes 
turning into a redshift, and the bisectors become more \-shaped. 
Each line shows its own individual behavior - function of its 
strength, excitation potential, and other properties. The line 
cores have been excluded from the bisectors for the Can lines. 
These wing bisectors show a progression from a C-shape to a 
\-shape from center to limb. 

Table |2] lists numerical values for the convective shifts com- 
puted as the mean of the bisector - excluding the part closest 
to the continuum - for each line and fi value, and Table [3] lists 
the line core velocities found when convolving the synthetic line 
profiles with the appropriate CRISP instrumental profile. 

4. 1 . Effect of spectral resolution 

The asymmetry of the lines will decrease with decreased spectral 
resolution. When comparing with observations, or using the the- 
oretical profiles for velocity calibration, the profiles must be con- 
volved to the proper spectral resolution. Figure |3] illustrates the 
dependence of the mean bisector, the center-of-gravity (COG) 
velocity, and the line-core velocity on spectral resolution. The 
core velocity is naturally the one most sensitive to spectral reso- 
lution. COG is least sensitive in this specific case, but the devia- 
tions in all cases start at such spectral resolutions that we are con- 
tent to use the mean bisector velocity in all our illustrations here. 
We do think that any practical use of our line profiles should use 
the detailed profile with the application of the appropriate instru- 
mental smearing and weighting. 

4.2. The limb effect 

The bisectors change with heliocentric angle, and close to the 
limb the blueshift usually gets small and is somet imes turned 
into a redshift for the line core s . This limb effec t ( lAdam et al.l 
II976I: iBeckers & NelsonI II978I: iBalthasad Il985h was subject 
to much spec ulation before granulation was understood (see 
lDravinsiri982|) . It seems that a comprehensive explanation is still 
lacking, which could be because the complicated nature of the 
situation precludes simple explanations. Nevertheless, we want 
to discuss it here in the light of our results. 

The convective blueshift and line asymmetry at disk center 
are caused by the asymmetry of large bright granules with a rel- 
atively slow upward flow producing blueshifted lines and nar- 
row dark intergranular lanes with fast downward flows produc- 
ing redshifted lines. As we leave the disk center, a number of 
effects will change the shifts and the line profiles. Towards the 
limb, higher photospheric layers are sampled that have different 
velocity patterns and different statistical relations between inten- 
sity and line-of-sight velocities. Of particular interest is the pat- 
tern of inverse granulation (e.g.. lCheung et alll2007l) that char- 
acterizes the photospheric structure above a certain height. Also, 
horizontal velocities become more important as the viewing an- 
gle decreases and will soon dominate the contribution from ver- 
tical convective motions. There are also 3D effects due to the 
corrugated nature of the optical depth surface. We can consider 
the following mechanisms that can decrease the blueshift and 
ultimately lead to a redshift in the average line profiles: 
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Fig. 2. Synthetic profiles and the corresponding bisectors resulting from the calculations of Ci, Fei and Can lines. The gray scale 
indicates the variation in the heliocentric angle from ^ - 1.0 to /i = 0.3. The bisectors for the Can lines are not shown in the core 
part. 



1. The intensity-velocity correlation changes as higher and 
high er levels in the ph otosphere are sampled by the line 
(e.g. lJanssen & Cauzzil2006.:, Cheung et alj2007h . This is the 
same effect that causes the C-shape in disk-center bisectors. 

2. The redshifted part of the granule is seen against a brighter 
background than the blu eshifted part (Beckers & Nelson 
ll978i:lAsDlund et aLllIOOOl) . 

3. Near the hmb, surface corrugation causes redshi fted parts 
of gra nules to show larger areas to the observer dBalthasaj 
[1981 . 

4. Reversed granulation "hills" over intergranular lanes shield 
blueshifted parts of granules when observed off disk center. 

Figure|4]shows a cross section through a simulation snapshot 
illustrating line formation at // = 0.3. Inspection of the simula- 
tion snapshots in this way can give support to all the listed mech- 
anisms and also inspire invention of new ones. The effect of the 
hot wall over an intergranular lane (at x ~ 2200 km) is, for ex- 
ample, obvious. It blocks the view towards the blueshifted side 
of a granule. But a similar structure is also present over the cen- 



ter of a granule (x ~ 4000 km), seemingly blocking a redshifted 
side of a granule. One can thus argue that some of the mentioned 
mechanisms and others can produce a blueshift with decreasing 
Deciding which mechanisms will dominate requires a quanti- 
tative analysis. We made a first effort in that direction. 

Effects |2] [3] and |4] are true 3D effects that break the symme- 
try between redshifted limbside granular parts and blueshifted 
centerside parts. They should go away in a multi-ID treatment 
where each column in the simulation is treated as a plane-parallel 
atmosphere. The results of such a numerical experiments for two 
of our lines can be seen in Fig. |5] It is clear that for each line, 
the variation of bisectors with /i is smaller than in 3D. It is also 
monotonous and always in the sense that a profile becomes red- 
der at the core of the line with decreasing //. The close following 
of the successive bisectors to that of /i = 1 is an illustration 
of the Eddington-Barbier relations: the /i dependence is equiv- 
alent to the depth dependence. We take this to mean that much 
of the center-to-limb variation of spectral line shifts is due to 3D 
effects. The latter either shifts the bisectors to redder or bluer 
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Fig. 3. Left: Bisectors of the synthetic Fei 6301.5 A line as a 
function of the spectral resolution. Right: Velocities measured in 
different ways from the profiles on the left panel, as a function 
of spectral resolution. 

wavelengths, relative to the ID case. For example, the blueshift 
relative to disk center that is seen in many lines for moderate 
heliocentric angles must be due to a 3D mechanism. We hypoth- 
esise that this comes from the screening of narrow downdrafts 
from view as the Une of sight departs from the vertical. 

4.3. Error analysis 

4.3.1. Non-LTE effects in Fei lines 

Our radiative transfer calculations assume LTE. The errors in- 
troduced by this approximation will to some extent be reduced 
by fitting the line profiles using a n individual abundance value 
for ea ch line as a free parameter. IShchukina & Truiillo Buenol 
(1200 ih studied the impact of non-LTE effects on Fei and Fen 
lines using a 3D simulation (but a ID treatment of the non-LTE 
problem, thus 1 .5D) and reported larger departures from LTE in 
granules than in intergranular lanes. Based on this work, we can 
expect the following two non-LTE mechanisms to dominate 

1. Overionization: A line from a neutral species typically has 
lower line opacity than its LTE value because of overioniza- 
tion that is caused by the UV radiation field being stronger 
than the local Planck function. This causes the line to be 
formed deeper in the atmosphere and thus be weaker than 
in LTE. 

2. The non-LTE line source function 5 J, can be expressed in 
terms of the departure coefficients of the populatio ns from 
the upper {bj) and lower {bi) levels ( see lRuttenl2003h if stim- 
ulated emission is neglected, 

S[ ~ ^BAT) (1) 

where By{T) is the Planck function. Photon losses in the line 
itself lead to a decrease in bj and increase on b,, producing 
a line source function that is lower than the Planck function. 
This leads to a stronger line than in LTE. 

The first effect tends to be the most important for weak lines, 
while the second is found for strong lines because the Une opac- 
ity dominates the background opacity. In some cases, the non- 
LTE and LTE profiles are very similar because the two effects 
cancel out. However, the presence of larger LTE departures in 
the granules than in the intergranular lanes not only changes the 
strength of the line, but also the convective line shift. We tested 
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Fig. 4. Properties of the model along a slice of one of the snap- 
shots at yU = 0.3; temperature, absorption coefficient at the core 
of Fei 630L5 A, and line of sight velocity from top to bottom. 
The solid line that is overplotted over the three panels represents 
the layer where tic — 1 along the line of sight. Four rays are 
shown as solid lines above t^q - 1 and dashed below this height. 
The bottom panel shows the line-of-sight velocities along the di- 
rection shown by the rays. The horizontal dashed line represents 
the layer where T500 equals 1 at ju = 0.3. 



the sensitivity of the resulting line shift to such differential depar- 
tures from LTE. One of the 3D snapshots was divided in granules 
and intergranular lanes based on the sign of the vertical veloc- 
ities at T500 = 1 and synthetic spectra were computed for dif- 
ferent values of the iron abundance. The aim was to find the 
combinations of abundances in granules and intergranules that 
reproduce a spatially-averaged spectrum with the same equiv- 
alent width as the observed profile. We used a nonlinear least- 
squares fit, described in Markwardt (2009), in order to infer these 
pairs of values. Figure |6] illustrates the average shift of the line 
as a function of abundance difference between the granules and 
the intergranular lanes. The reference was set where the abun- 
dance derived from the granules and intergranules is the same. 
Shchukina & Trujillo Bueno (2001) find an average non-LTE ef- 
fect on derived Fe abundances of 0. 1 dex. If this number is taken 
as also typical of the difference between granule and intergran- 
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Table 2. Convective shifts in m s ' of the studied lines from ju - 1 .0 to /i = 0.3. 



Element 


^0 (A) 


yU = 1.0 


11 = 0.9 


11 = 0.8 


11 = 0.7 


= 0.6 


11 = 0.5 


11 = 0.4 


11 = 0.3 


Ci 


5380.34 


-948 


-936 


-922 


-910 


-897 


-888 


-880 


-844 


Fei 


5250.21 


-294 


-372 


-410 


-426 


-420 


-394 


-346 


-259 


Fei 


5250.65 


-304 


-354 


-370 


-369 


-348 


-308 


-245 


-140 


Fei 


5576.09 


-334 


-378 


-388 


-383 


-360 


-317 


-249 


-140 


Fei 


6082.71 


-441 


-495 


-522 


-536 


-531 


-511 


-473 


-397 


Fei 


6301.50 


-307 


-349 


-358 


-351 


-325 


-279 


-211 


-101 


Fei 


6302.49 


-335 


-387 


-404 


-403 


-381 


-342 


-280 


-176 


Fei 


7090.38 


-352 


-396 


-409 


-407 


-389 


-354 


-302 


-210 


Can 


8498.01 


-265 


-288 


-289 


-275 


-248 


-206 


-148 


-70 


Can 


8542.05 


-335 


-364 


-374 


-372 


-359 


-337 


-307 


-269 


Can 


8662.16 


-272 


-294 


-293 


-279 


-251 


-208 


-149 


-68 



Notes. The bisector of each line was integrated from 95% of the continuum intensity to the core of the line. In the case of the Can lines, the 
non-LTE core was excluded. For C i the upper threshold was higher (97%) because the line is weak. 

Table 3. Line core velocities of the synthetic spectral lines (excluding Ca ii) convolved with the appropriate theoretical profile of 
CRISP for each wavelength. 



Element 


^ (A) 


ji = 1.0 


11 = 0.9 


11 = 0.8 


11 = 0.7 


11 = 0.6 


fi = 0.5 


11 = 0.4 


fi = 0.3 


CI 


5380.34 


-928 


-895 


-867 


-840 


-820 


-810 


-795 


-764 


Eel 


5250.21 


-178 


-306 


-358 


-370 


-354 


-315 


-246 


-134 


Eel 


5250.65 


-115 


-196 


-221 


-213 


-177 


-115 


-26 


108 


Eel 


5576.09 


-163 


-222 


-231 


-209 


-164 


-96 


-3 


126 


Eel 


6082.71 


-403 


-470 


-499 


-508 


-500 


-480 


-433 


-341 


Eel 


6301.50 


-188 


-238 


-241 


-217 


-173 


-109 


-17 


108 


Eel 


6302.49 


-255 


-326 


-347 


-338 


-306 


-255 


-178 


-57 


Eel 


7090.38 


-315 


-378 


-396 


-388 


-362 


-321 


-255 


-148 




Velocit}^ [m • s ^] Velocity [m • s ^] 



Fig. 5. Bisectors for two lines computed in the plane-parallel 
treatment, "ID", compared to the full "3D" result. 



ule abundance in Fig.|6] the errors in the computed lineshift due 
to non-LTE effects should be smaller than +30 m s"' . 

4.3.2. Profile fitting 

Due to the effects described in Section l4.3.1l and becuase the 3D 
models are only realistic up to a certain height, the resulting pro- 
files show differences compared to the FTS atlas. The solar abun- 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

A(granules)-A(inter-granules) [dex] 



Fig. 6. Differential Doppler shift of the Fei 6082.7 A line pro- 
duced by changes in the relative abundance between granules 
and inter-granules. The combinations of granular/intergranular 
abundances were fitted to reproduce the equivalenth width of the 
FTS atlas. 

dances inferred bv lAsplund et al.l (l2005h normally lead to differ- 
ent results for the different lines, some presenting more damped 
wings than the atlas and others narrower ones. The abundance is 
used as a free parameter in order to match the wings of the syn- 
thetic profiles for ii - I with the FTS atlas. The inner core is not 
included in the fitting procedure becaus e generally non-LTE ef - 
fects are expected to be larger there (see lRutten & Kosti5ll982h . 
and the model snapshots become le ss reliable close to the upper 
boundary (see lAsplund et alJl2000l) . 

As an indicator of the errors due to this fitting procedure, we 
performed test calculations to get the sensitivity of the bisectors 
to changes in the assumed abundance. The results are shown in 
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Fe I 5250.2 A Fe I 5250.6 A Fe I 5576.1 A Fe I 6082.7 A Fe I 6301.5 A Fe I 6302.5 A Fe I 7090.4 A 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

-600-400-200 -600-400-200 -600-400-200 -600-400-200 -600-400-200 -600-400-200 -600-400-200 
Velocity [m • s^'] Velocity [m ■ s^^] Velocity [m ■ s^'] Velocity [m ■ s^^] Velocity [m • s^'] Velocity [m ■ s^^] Velocity [m ■ s^^] 

Fig. 8. Fe I synthetic bisectors (fuUdrawn) compared with bisectors measured from the FTS atlas (dashed). 




-300 300 600 -1060 -900 -750 -600 

Velocity [m ■ s"-*-] Velocity [m • s"-*-] 




A(Fe) [dcx] A(C) [dex] 



Fig. 7. Bisectors calculated for at = 1.0,0.6,0.3 for a range 
of abundances. The results are only from five snapshots which 
means that they cannot be directly compared with, e.g., Fig|2] 
Left: abundance test for Fe i 6301 A with values of the Fe abun- 
dance within [7.25-7.65] . Right: test for Ci 5380 A, the values 
of the abundance span an unrealistic range compared to the nom- 
inal value [8.19-8.69]. 



Fig- 13 The left panels correspond to the Fe i 6301 A line, calcu- 
lated on three different heliocentric angles, for different values 
of the Fe abundance. The panels on the right correspond to the 
C I 5380 A, where the abundance covers an unrealistic range of 
values [8.19, 8.69]. This weak line, which is formed very close 
to the continuum, is more sensitive to changes in the abundance 
than the stronger Fe i line. The lower panels show the relative ve- 
locity shift integrated over the whole bisector We estimate that 
the uncertainties from the fitting procedure are within 0.1 dex 
in abundance which then corresponds to +10 m s ' for the inte- 
grated bisector. 

4.3.3. Discussion 

From the tests above, the estimated error in the shift of the lines 
is less than 50 m s This means that a full non-LTE treatment 
combined with the exact atomic data should not change these re- 



sults by more than this amount. The estimate does not, however, 
include the intrinsic errors of the simulation. 

To check the realism of the simulations, we compared the 
synthetic Une profiles at = 1.0 with the FTS atlas after ap- 
plying a gravitational redshift on the former. The average syn- 
thetic Ci bisector lies 530 m s"' redwards of that measured 
from the FTS atlas, something that we deem mostly to be due 
to error in the laboratory wavelength of the line. At face value, 
our simulations implies a wavelength of 5380.3275 A on the 
FTS-atlas wavelength scale. This is within the cited uncertainty 
0.02 A ) of the laboratory value of 5380.3370 A from lJohanss"onl 
T966I) . which was used here. It is also comfortingly close to the 
5380.3262 A found in a similar way bv lLangangen et al] (l2007h 
from their observations and simulations of magnetic regions. 
The bisectors for the Fe i lines are displayed in Fig. [8] The syn- 
thetic bisector shapes show an acceptable correspondence with 
those measured from the FTS atlas. The standard deviation of 
the differences in mean velocity is 56 m s"' which corresponds 
to approximately 1 mA. This is comparable to our estimate of the 
internal errors from non-LTE eff'ects and fitting procedures, but 
once again there is strong reason to believe that a significant part 
of these deviations are due to errors in the adopted laboratory 
wavelengths. 

Finally we note that lAUende Prieto & Garcia^ Lopezl(ll998h 
estimate the precision in FTS wavelengths to be on the order 
of 50 m s"'. It seems likely that several of the error estimates 
discussed above are too conservative, but the synthetic profiles 
seem to pass the test and they should be correct within 50 m s"' . 
Observational checks of the off-center profiles are more difficult 
but would be welcome. Since photospheric modeling often be- 
comes more uncertain for low ji values, it is natural to suspect 
that the accuracy of our convective line shifts will decrease to- 
wards the limb. 

The hydrodynamic granulation simulations do not include 
magnetic fields. This is potentially the most important source 
of errors in the practical application of these s ynthetic line pro- 
files f or velocity calibrations. For example. 'Langange n et al.l 
(l2007l) estimated a difference of 200 m s ' in convective shift 
for the C i 5380 A line between non-magnetic and active gran- 
ulation. (This was mostly caused by the changing intensity and 
velocity pattern and not by Zeeman splitting.) Any calibration 
work must use the least magnetic photospheric region avail- 
able. However, how strong the quiet Sun ma gnetic fields are 
is still an ongoing debate (e.g. [M artinez Gon zalez et al J [20081 : 
iPietarila Graham et alJ[2009[:[deWiin et al.li2009lK 
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4.4. Calibration of SST/CRISP observations 

A quiet Sun dataset acquired with SST/CRISP at = 1.0 
was calibrated using the current results. The Fe i 6302 A line 
was sampled with eight wavelength points, evenly distributed 
with a separation of dA = 48 m A and processed usin g the 
image restoration code MOMFBD dvan Noort et al. I |2005l) . The 
telecentric optical design of CRISP (Scharmer l2006ir IS opti- 
mum for high image quality, but it introduces instrumental field- 
dependent wavelength shifts produced by irregularities on the 
surface of the etalons, which are commonly known as cavity er- 
rors. The latter complicate the task of obtaining the spatially- 
averaged spectrum, as the profile is effectively broadened by 
these random pixel-to-pixel shifts. The average spectrum was 
computed with a least-squares-fit of all the data points to a spline 
which was characterized with ten parameters. For this purpose, 
all the individual spectra were corrected for cavity errors and 
then arranged in a large array, which was sorted in increasing 
wavelength. 

To carry out the velocity calibration, the synthetic spectrum 
obtained from the 3D simulation was convolved with the the- 
oretical instrumental profile of CRISP, which included the ef- 
fects of the beam converging at //165. The velocity calibra- 
tion was determined by finding the constant wavelength shift 
to the observed bisector that minimized the difference to the 
synthetic one. Figure |9] illustrates the results of the calibration, 
where the Doppler map at the core of the line is scaled between 
+ 1500 m s The bottom-left panel contains a density map of 
the spatially resolved spectra, corrected for cavity errors. The 
dispersion of points around the seven nominal wavelengths ap- 
pears after correcting for the cavity errors in combination with 
the granulation contrast. The bisector of our calibrated data is 
overplotted to the bisector of the synthetic spectrum that was 
convolved with the instrumental profile of CRISP. There is a tel- 
luric blend on the red wing of the observed profile, not present 
in the synthetic spectrum, that pollutes the upper part of the bi- 
sector. 



5. Usage of electronic data 

The synthetic profiles that have been presented in this study can 
be accessed electronically. The spectral response of any instru- 
ment produces broadening and can introduce asymmetries in the 
observed profiles. Thus it is imperative that the profiles be de- 
graded with the proper instrumental response before they are 
used to calibrate observations. 

The quiet Sun's average spectrum can be computed from 
quiet regions in the observations. However, the presence of p- 
modes (oscillations) will produce an unrealistic shift in the ob- 
served average profile. Thus, it should be averaged in time and 
space, or at least from a sufficiently large sample of quiet Sun 
spectra. 

We have been cautious in the treatment of line cores, but the 
comparison with the FTS atlas suggests that the errors in the 
cores are small. Except for the Can IR triplet, we do not dis- 
courage including the cores in velocity calibrations. We think, 
however, that a velocity measure that weigh in more of the pro- 
file will give a more reliable result. 

6. Conclusions 




-0.12 -0.06 0.00 0.06 

Wavelength [A] 



Velocity [m • s 



We m ade use of the approach described by (iLangangen et alJ 
I2OO7I) to obtain an absolute reference for velocities in solar ob- 
servations. For this purpose, we calibrated the shift in synthetic 



Fig. 9. Top: Doppler map of quiet sun granulation observed with 
SST/CRISP at Fe i 6302 A. The image gray scale is clipped to 
+ 1500 m s"'. Bottom-left: Spatially-averaged profile computed 
along the image shown in the top panel. The density map rep- 
resents the spatially-resolved profiles corrected for instrumental 
cavity errors. The average profile was computed using a least 
squares fit of all the points to a spline, consisting of 10 param- 
eters. Bottom-right: Bisector of the observed spatially-averaged 
profile. The velocity calibration was carried out using our syn- 
thetic profile convolved with the CRISP instrumental profile. 



spectra produced by convective motions at different heliocentric 
angles for C i, Fe i, and Ca 11 lines as commonly used for solar 
physics diagnostics. We propose to use the profiles as absolute 
local references for line-of-sight velocities to allow accurate ve- 
locity measurements. 

A realistic 3D hydrodynamical simulation of the non- 
magnetic solar photosphere was used to synthesize the spec- 
tra. The computations were done in LTE for C i and Fe i lines 
whereas non-LTE computations were carried out in the Ca 11 in- 
frared lines. As the conversion from wavelength to velocity scale 
is done using the same laboratory wavelength as used for com- 
puting the profiles, the measured line shifts are relatively unaf- 
fected by errors in the central wavelength. Uncertainties in the 
rest of atomic parameters are partially compensated when the 
profiles are fitted to reproduce the strength of the FTS atlas. 



8 



J. de la Cruz Rodriguez et al.: Solar velocity references from 3D HD photospheric models 



Our results show a systematic center-to-limb variation of the 
bisectors, which is determined to a large extent by the 3D struc- 
ture of the atmosphere, which is projected along the line-of- 
sight. However, there is a small amount of redshift (relative to 
jU = 1 .0) produced by the change in the formation height with 
the heliocentric angle (non-3D effects), which is not related to 
projection effects. A careful analysis of the models suggests that 
a large number of 3D snapshots must be used in studies of line 
shifts so the oscillations are properly sampled and the statistics 
of granules and intergranules are complete, whereas non-LTE ef- 
fects seem to have a minor influence on the shift of our Fe i lines. 
The largest uncertainty is related to the non-magnetic nature of 
the simulation. Nevertheless, we estimate our calibrations to be 
accurate up to ~ 50 m s"' at = 1, remaining untested towards 
the solar limb. 

We would like to emphasize the somewhat limited usability 
of the Can calibration, as observations are polluted by a sig- 
nificant number of blends, which should be masked in order to 
compare with the calibration data. 

Our synthetic profiles reproduce spatially-averaged solar ob- 
servations accurately, but it would be interesting to test for the 
impact of magnetic fields on the convective shift of the lines, and 
confirm the reliability of our calculations near the limb. 
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